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SINGULAR LIMITS IN PHASE DYNAMICS 
WITH PHYSICAL VISCOSITY AND CAPILLARITY 

K.T. JOSEPH AND P.G. LEFLOCH 

Abstract. Following pioneering work by Fan and Slemrod who 
studied the effect of artificial viscosity terms, we consider the sys- 
tem of conservation laws arising in liquid-vapor phase dynamics 
with physical viscosity and capillarity effects taken into account. 
Following Dafermos we consider self-similar solutions to the Rie- 
mann problem and establish uniform total variation bounds, al- 
lowing us to deduce new existence results. Our analysis cover 
both the hyperbolic and the hyperbolic-elliptic regimes and apply 
to arbitrarily large Riemann data. 

The proofs rely on a new technique of reduction to two coupled 
scalar equations associated with the two wave fans of the system. 
Strong convergence to a weak solution of bounded variation 
is established in the hyperbolic regime, while in the hyperbolic- 
elliptic regime a stationary singularity near the axis separating 
the two wave fans, or more generally an almost-stationary oscil- 
lating wave pattern (of thickness depending upon the capillarity- 
viscosity ratio) are observed which prevent the solution to have 
globally bounded variation. 



1. Introduction 

The Navier-Stokes equations for van der Waals fluids with viscos- 
ity and capillarity effects included, allow one to model the dynamics 
of liquid-vapor flows. The associated set of first-order conserva- 
tion laws is of hyperbolic or hyperbolic-elliptic type, and admits 
propagating discontinuities (shock waves). The singular limit cor- 
responding to vanishing viscosity and capillarity coefficients allows 
one to select physically admissible, discontinuous solutions to the 
first-order conservation laws. In particular, this yields an approach 
to select solutions to the so-called Riemann problem, when the ini- 
tial data consists of two constant states separated by a single jump. 
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Our purpose in the present paper is to derive uniform (viscosity- 
capillarity independent) bounds and to justify this singular limit in 
the context of nonlinear elasticity and phase transition dynamics. 

Attention will be concentrated on the following system of two 
conservation laws: 

— - —(o(w) - e — + 5^^] - 

dt dx^ dx dx^ ' ' ^YY) 

dw dv _ ^ 
dt dx 

where v = v{t, x) and w = iv{t, x) represent the velocity and deforma- 
tion gradient of the fluid or solid material under consideration. We 
consider also the associated first-order system 



dW _ g 

dt dx 



(1.2) 



and impose Riemann initial data 



{v,w){0,x) = < (1.3) 

[Vr, IVr, X> 0, 

where Vi, Vy, Wj, w ^ 6 R are given constants. Solutions of (|1.2 |) - (|1.3|l are 
known to be self-similar, that is, to depend only upon the variable 
y := x/t. 

When viscosity and capillarity terms are taken into account, the 
corresponding set of differential equations reads 

-yv' -o{wy = ev" -6w"', 

(14) 

-yw'-v' = 0, ^ ^ 

supplemented with the boundary conditions 

lira {v,w){y) = {vi,ivi), \im {v,w){y) = {Vr,iVr). (1.5) 

I/— >— CO y— >+oo 

Our main objective in this paper is to establish : 

(1) the existence of a smooth, self-similar solution to (|1.4 |) - (|1.5|) 
having uniformly bounded total variation, 

TV{Ve, We) < \Vr - Vl\ + \Wr - Wl\, (1.6) 

with the implied constant being independent of e, 6 within the 
range << 1, and 

(2) the strong convergence of Ve, We toward a weak, self-similar 
solution V, w to the first-order conservation laws (|1.2|) . 
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An outline of our main results is as follows. 

Assuming first that the system (|1.2[) is hyperbolic we provide a 
rather direct proof of the above two properties, first in the viscosity- 
only case (Section 2) and then for general viscosity-capillarity (Sec- 
tion 3). In this regime, it is known that the limiting solutions in 
general contains nonclassical shocks, with depend upon the ratio 
see d. 

Next, we investigate the generalization to the hyperbolic-elliptic 
regime, and discover a concentration phenomena near the axis x = 
separating the two wave fans. This seems to be consistent with 
numerical experiments with the model under consideration, but it 
would be interesting to check numerically the feature discovered here 
analytically. In the context of phase dynamics, it is also well-known 
that the limiting solutions contains subsonic phase boundaries which 
again depend on the capillarity to viscosity ratio; see l|23l |25l l27l [Ij 

Observe that our results cover arbitrary large Riemann data and 
physical viscosity and capillarity terms. Our results supplement the 
earlier, pioneering work by Fan and Slemrod [7, 8J where an artificial, 
"full" viscosity was used. Our technique of proof in the present paper 
is quite different from the one in [7J, as we introduce a decomposition 
of the system of equations (|1.4|) into two coupled scalar equations. 

Finally, we provide various generalization to the boundary value 
problem and more general classes of conservation laws. 

Recall that the activity on self-similar vanishing viscosity limits 
started with an extensive research by Dafermos t2l-|5l (see also [6J) 
who advocated the use of self-similar regularizations to capture the 
whole structure of wave fans within solutions of the Riemann prob- 
lem. Self-similar approximations in the context of phase transition 
dynamics were studied by Slemrod ||23ll25l . Fan and Slemrod [7], who 
covered large data solution and artificial regularization terms. Next, 
small data solutions to general systems were treated by Tzavaras ||28| 
(conservative systems) and LeFloch and Tzavaras |(T8i [191 (noncon- 
servative systems). 

Self-similar diffusive-dispersive approximations for general sys- 
tems were investigated in LeFloch and Rohde fV7\. Joseph and 
LeFloch |9l-|T2l extended the technique to cover boundary value 
problems, relaxation approximations, and general diffusion matri- 
ces. The present paper is the continuation of Il9l-|[T2|. For issues 
related to the discretization of the viscosity and capillarity terms, we 
refer to BMIll. 
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2. Elastodynamics with physical viscosity 

We begin with the case that 6 = and Oi„ > 0. We will prove: 

Theorem 2.1 (Vanishing viscosity limit in elastodynamics). Assume 
that the system (|1.2|l is uniformly hyperbolic, in the sense that there exits a 
positive constant Cq such that 

inf o^v{iv) > Cq > 0. 

iveR 

Then, given arbitrary Riemann data vi, Wi and v,-, Wr, the viscous Riemann 
problem (|1.4 [) - (|1.5|) with 6 = admits a solution v^, which has uniformly 
bounded total variation and converges pointwise to a limit v, w. This limit 
is a weak solution of the Riemann problem (|1.2|) - (|1.3|) . 

It is worth observing that the solutions w^, v^ contain a mild singu- 
larity at y = 0: their derivatives up to order C/e only are continuous 
at y = 0. Hence, this singularity vanishes in the limit e — > 0. The 
singularity is due to a factor y appearing in the key equation (|2.2|) 
below, in front of the term with the highest order derivative. 

Proof. Step 1. Reduction to a scalar equation for the component w. To 
simplify the notation we suppress the subscript e. We will first study 
the problem on a finite interval [-L, L] with L chosen to be sufficiently 
large and we will show the existence of a solution to (|1.4[) satisfying 
the boundary conditions 

{V,w)i-L) = {Vi,Wi), {V,W){L) = {Vr,Wr). (2.1) 

First, we observe that the system can be reduced to a single scalar 
equation for the unknown w. Namely, using v' = -y w' in the first 



equation of (|1.4|) we find 

(y^ + e - o,a{w)) w' = -ey w" , (2.2) 

which will be studied in two regions [-L, 0] and [0, L] with the fol- 
lowing boundary conditions 

w{-L) = wi, w(O-) = w», (2.3) 

w(0+) = w„ w{L) = Wr, (2.4) 

where w^ is a parameter to be determined so that the boundary con- 
dition on the variable v (which has been eliminated from (|2.2[) ) is 
satisfied. This condition is obtained by integrating over [-L, L] the 



second equation in (|1.4|) : 



"Vr - ^/ = I v' dy = -Lwy + Lwi + | w dy. (2.5) 

J-L J-L 
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This is a global, integral condition on the solution w. We will now 
solve the problem (|Z^ -(|23]). 

Our boundary condition can be justified as follows. If we search 
for a piecewise smooth solution of (|1.4|) , then clearly the measure 
v' = -yiv' can not have a mass point at y = 0, hence v is continuous. 
Similarly, by the first equation in (|1.4[) the function o{w) — ev' must be 
continuous at y = 0. But, using again the equation ev' = -yew' we 
see that a{w) itself must be continuous. Finally, since a is assumed to 
be strictly increasing, w must be continuous. 

Based on the function w we can define a{y) := y - (ai„(w(y)) - e)/y, 

^a{x)dx on [-L,0) and on (0,L] with a~ e 
[-L, 0) and a'^ G (0, L], respectively. Let us show that the minimum 
with respect to y of this integral is attained at 

:= ±{aUw{p^))-ef'. 
For definiteness, consider the case a = a~ e [-L,0) with a~ < y; then 

a{x)dx> {x-^-^ -)dx>{y^ -a^)l2+\o^{alyy^-', 

J a J a 

which tends to infinity as y — > 0. Necessarily, the minimum is at- 
tained away from 0. As this integral has quadratic growth in y, for 
L large the minimum p = p" is attained in (-L, 0) and is given by the 
equation p^ - Oiu{w{p)) + e = 0. 

This choice of p- gives a{x) dx > 0. By setting 



/pi a{x) dx -{lie) f^^ a{x) dx 



the problem under consideration is equivalent to 



w{y) 



wi + (w, - wi) I (p- dx, y < 0, 

^± (2.6) 

Wr + (W, - Wr) I (p+ dx, y > 0, 

"^y 

together with the boundary condition (|2.5|) . 

To find w„ we used the continuity of the solution at y = namely, 
w{0-) = w{0+) - w* and v{0-) = v{0+) = u*. Integrating the equation 
v' = -yw' from -L to and from to L and using (|2.6|l , we get 



- = (w» - wi) I -y(p-{y)dy. 
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Vr-v, = {tv, - Wr) \ x/(p+{y)dy. 
Jo 

Adding these formulas and solving for zv*, we get 

Vr-vi + Wr y(p+dy - Wi y (p- dy 



f^yq^idy- £^y(ptdy 



(2.7) 



which provides us with the value of w». 

Observe that the denominator is bounded away from zero. To 
evaluate its minimum value (in terms of the constitutive function o 
and the Riemann data) we now need to analyze the functions (p+. 

Step 2. Properties of the functions (p+. We will show now that the 
support of the functions <p+ is essentially concentrated away from 
the axis y = 0. Let o"^ = min_L<y<L o^{w{y)), = max_L<y<L Ozu{w{y)), 
= {oZ - e)i/2 and Af = {o^ - eY'^. We claim that there exists a 
constant C > such that 



< (p_(y) < - \ 



2e 



[{-lylX^)- 



-L<y< -Af , 
-Af<y< -Ar, 
-Af < y < -Af /4, 
-A-/4 < y < 0, 



(2.8) 



< (p+(y) < - 



(2y/An^, 0<y<Ar/4, 
Af /4 < y < Af, 
A- < y < Af, 
Af < y < L. 



2e 



2e 



(2.9) 



This is proved as follows. First, in the interval < y < A^/4 we 
find 



Jo+ ^ Ja™/2 ^ 



dx > I dx 

A?/2 

Jy X 



/9i/^(3/4)A,^ 



log(A-/2y) 
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and, next, for y > 

(x^ - o-a,ix) + e)dx > — -dx 

Jp+ Jam X 

= Cix - Mx) + ef'){l + ^^"^"^^ ^ ^^'" )dx 
Jam X 



L 



y (y - Aff 

X — )UX = 



> (x-Af)dx = 



The other cases are completely similar. Since, on the other hand, it is 
easy to check that 

this completes the derivation of the properties of (p+. 

Step 3. Existence result. Observe that w, consists of the sum of 

Vr -V; Vy- Vi 



Jo y <p+ - u y "p- Lf -y ^y + U yp^^y^ 

and a convex combination of iVi,iVr. In view of this observation, 
together with the estimates (|2.8|) and (|2.9|) , we get 

\Vr - Vi\ 



\w^\ < max(w/, w,.) + 



2Af + 0(e") 



\Vr - Vl\ 

< max{wi,Wr) H =: Aq 

Co 

and we arrive at the uniform bound 

Kl < Aq. (2.10) 

Replacing now (by its value given in (|2.7|l ) in the formula (|2.6|) 
we can define a mapping w 6 C°([-L,L]) i-^ r(w) G C°([-L,L]). We 
claim that, for fixed e, the function r(w) is of class and 

||r(a;)||co < Ao, 
C 

Il'r'(w)'llc0 ^ ~ (l^r ~ W»l + 1^1- ~ 

e 

Indeed, these estimates follow immmediatly from (|2.6|) , (|2.8|) , and 

In consequence, the operator T is a compact map from the convex 
bounded set {w / ||w||co([-l,l] ^ Aoj into itself. By Schauder's fixed 
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point theorem, T has a fixed point which satisfies (|2.6|) and, clearly, 
is of class C^. Furthermore, in view of (|2.6)) and (|2.7[), we have 



I 



2 

W'{y)\dy < \Wr - W*| + |W/ - W*| < |Wr - W/l H \Vr - Vl\ (2.11) 

L Co 



and 



J b'(y)|dy< J \yw'{y)\dy <{Af + V)i^Wr-wi\ + j^\Vr-vi^. (2.12) 

Next, since all of the estimates are uniform in L, we can let L — > oo 
and we conclude with existence of a solution {v'^ , w^) on defined on the 
whole real line (-oo, oo). This function (u'^, w*^) is a bounded solution 
of (|1.4[) - (|1.5|) and has uniformly bounded total variation 



%J —a 



{^w'{y)\ + \v'{y)\) dy < {Af + 2) (\w, - + ^ \Vr - Vi\). 



By Helly's compactness theorem, it admits a subsequence converging 
pointwise to a limit {v, w), as e goes to zero. Clearly, the limit is a 
weak solution of the problem (|1.2 |) - (|1.3|) . □ 



3. Elastodynamics with physical viscosity and capillarity 

In this section, we consider the full system with physical viscosity 
and capillarity included. We will prove: 

Theorem 3.1 (Vanishing viscosity-capillarity limit in elastodynam- 
ics). Assume that mfo-,o > Cq > 0. Then, given arbitrary Riemann data 
Vi, Vr, lUi, lUr, the problem (|1.4|) - (|1.5|) with 6 = ye^, y > 0, admits a solution 
v^, which has uniformly bounded total variation and converges to a limit 
V, w. This limit is a weak solution to the Riemann problem (|1.2 [) - (|1.3|) . 

Observe that this is a large data result. 

Proof. Step 1: reduction to a scalar equation. We will first study the 
problem on a finite interval [-L, L] with L chosen to be sufficiently 
large with boundary conditions 

{V,W){-L) = (VuW,), {V,W){L) = {Vr,Wr). (3.1) 

As before, the system can be reduced to a single scalar equation 
for the unknown w. Namely, using v' = -y w' in the first equation of 
(iLill we find 

{y^ + e- a^{w)) w' = -ey w" - bw'" , (3.2) 
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which will be studied in two regions [-L, 0] and [0, L] with the fol- 
lowing boundary conditions 

w{-L) = Wi, w{0-) = w„ (3.3) 

w{0+) = w„ w{L) = Wr, (3.4) 

where is a parameter to be determined so that the boundary con- 
dition on the variable v (which has been eliminated from (|2.2[) ) is 
satisfied. This condition is obtained by integrating over [-L,L] the 
second equation in (|1.4|) : 

Vr -Vi = J" v' dy = -L zvr + Livi + w dy. (3.5) 

This is a global, integral condition on the solution w. We will now 
solve the problem (|3.2|) - (|3.5|) . 

We will follow a method originally introduced by LeFloch and 
Rohde I fTTI for small amplitude solutions of systems of conservation 
laws. Setting (p = e^H, the differential equation under consideration 
becomes 

6[H" + 2H'|S' + H|S" + H|S'2] + ey[H' + H^'] + {y^ + e - o,,)H = 

Setting the coefficient of H' equal to 0, we get jS' = thus jS(i/) = 
and 

H" = tH^a^^H. (3.6) 



where 

,2 



2 

e _ e 



lj{zv{y), y) = Ou, + y\— - 1) - -. (3.7) 
We take 6 = ye^ and, in this case, 

1 e 
li{w,y) = Ou, + y\— - 1) - - > 0, 

provided e, y > are sufficiently small. Upon integrating this equa- 
tion for H and substituting, we get 



where 



li-^l\fdx+ ii-^l\"\dx 

L J-L 
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and 



. r^z^ + (3.9) 

Jp 2y y y 
= ^ I + Jl + ^y{Oza - e/2 - x^))dx. 

Here, p is the maximizer of the the integral 



r\-x hi{wix),x) 
max I (- — h )dx. 

-L<y<LX 2y A/ y 

Indeed p satisfies 

p^ - Ora{w{p)) + e/2 = 0, 
which has pairs of solutions p_ G [-L, 0) and on p+ e (0, L], satisfying 

We consider cp- and (p+ on [-L, 0] and [0, L], defined by 

1 + 0(y) iKy^ 

^^^^^^ = (i^^ ^ ^'-^'^ 



and setting 



, , ^-(y) , , ^+(y) 



j_^(p-{x)dx (p+{x)dx 
the problem under consideration is equivalent to 

{wi + (w, - wi) (p- dx, y <0, 
\ . . (3-12) 

+ (^r - w,) (p+ dx, y > 0, 

together with the boundary condition (|2.5|) . The latter can be ex- 
pressed in its original form involving w' (see the second equation in 
((L4|) ) and yields 

Vr-vi + Wr C w+^y ~ ^/ y *P- ^y 

:= 1 5 . (3.13) 

X y(p+dy- l^ycp-dy 
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Observe that the denominator in the formula for iv* is bounded away 
from zero. To find the value of the minimum in terms of the con- 
stitutive function o and the Riemann data we need the following 
properties of (p± 

Step 2: Properties of (p+. With the same of and as before, set 

Af = {oZ-e/iy^^andA'^ = {o^-e/iy^. There exist constants Q > 0, 
and C > depending only on y, o^, and such that 



< (p_(y) < ^ 



e , 
1, 

^ € ' 



-L < y < -\f, 
-Af < y < -Af , 
-A^ < y < -yi/^A^, 

_yl/2^m < y < 0, 



(3.14) 



< <p+(y) < 



C(y-A™)2 



- C(y-Aff 



< y < 7^/2 

K<y< Af , 

A^<y<L. 



(3.15) 



First, we prove the estimates for (p+. Consider the region y > A^ 



, , r (-X L-e/2 + x2(l/4y-l). 

IT/-, /^i , (cT^-e/2-x2), 



27 Ja*^^ 



1+ 



(cr„; -ell- x^) 



^dx, 
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and, for y > , 

Jam 



Af r 



2 



<— (Af-x)ix = - 



y jAf y 2 

Now, consider the region 7^^^A^/2 < y < A^ , 

, , r i-x k-e/2 + x2(l/4y-l). 
^^^'^^^ ^ Jp. W ^ \ 7 ^ 



— (i-ji+47 — - — yx, 



and, for y^^^A"J/2 < y < Af , 

p(y,p.)-<^J^ 47 )dx 

' 47(c7^ - e/2)i/2 _ _ g/2)i/2 + x)dx 

y 

<-C j {a^-e/2f'^ -x)dx 
Jy 

<-C I (A^-:c)dx = -C^^^-— 



-J 
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For < y < y^'^A'^^ll, 

, . r i-x | g.„-e/2 + x^(l/4y-l) ^^ 



(-_x \o,,-m^x^my^.^^ 



2y yi/2 



= I -^*=-ciy-'>"')'"'i- 

The other cases are identical. One can also check (for some c > 0) 
I (p.{y)dy>ce, I (p+{y)dy>ce, 

J-L Jo 

SO that the desired properties of (p+, cp-, follow. 

Existence arguments. Now w, is sum of convex combination of 
ivi, Wr and the quantily 

Vr-Vl _ Vr - Vl 

^ycp^y- U y "P- Lf -y(p-dy + ^ y(p+dy + 0{e-) 

From this we get 

< max(w,,w,) + ^^„, ^ < Ao. 

We can now replace (given by (|3.13|) ) in \3.12) and arrive at a 
mapping w G C°([-L,L]) T(w) g C°([-L,L]). For fixed e, T{w) is of 
class C^, and 

||r(K;)||co < Ao, 

c 

l|T(w)'||cO < — {\Wr - W,| + \Wr - 

e 

These estimates follow from (|3.12|) (|3.14|) and (|3.15|) 

Thus, T is a compact map from a convex bounded set {w : \\iv\\co < 
Aq} into itself. By Schauder's fixed point theorem, T has a fixed 
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point in C°[-L,L] and satisfies (|3.12|) and so the solution is of class 
C^. Furthermore, in view of (|3.12[) and (|3.13|) , 

w'{y)\dy < \Wr - w*| + |w/ - w,| < |Wr - W/l + 1}— — ^ (3.16) 

L Co 



x: 



and 



\v'{y)\dy < J |yw'(y)|dy < (A^ + 1) - wi\ + t"' ^ ^'' ). (3.17) 

As the estimates are uniform in L, we have existence of solution 
(^^ ,w^) on (-00,00) we have (;p^ ,w^) bounded solution of (|1.4|) and 
(|1.5|) with uniform total variation 



r 

%j —o> 



{\w'{y)\ + \v'{y)\dy < (A^ + 2)[\Wr - wi\ + 2^^^-^] 

Co 



By compactness there exists a sequence converges in and pointwise 
a.e. to a function (v, w) as e goes to zero and solves (|1.2[) and (|1.3|) 
with e = 0. □ 

4. Phase dynamics with physical viscosity 

We turn to the model of phase transition dynamics when the system 
is not strictly hyperbolic, that is, when Ozu{w) is only non-negative, or 
even is hyperbolic-elliptic when Ojg takes negative values. 

To illustrate the key difficulty we will have to cope with, let us first 
consider the example 

4y) = y + c/y, 

log(|y/p+n and so 
In this case 

e^^^iylb)-'!' 



(p+(y) 



^\~-^{ylb)-^l^dy 
which is concentrated at y = 6. For c < 0, we have, p+ = (-c)^/^ and 

-if a{x)dx dilll^. sXllx-cle 

SO that 

e^(y/(-c)l/2)-c/^ 



/;^e^(y/(-c)V2)-./.dy 
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which is concentrated at y = (-c)^^^. 

In the previous section, to show that is bounded independent 
of e it was crucial that (p* and (p~ be concentrated away from 0. In the 
above example, if c > 0, we have seen that this is not the case. This is 
the difficulty is getting estimates when au,{w{y)) oscillates, i.e. takes 
both negative and positive values. 

We now prove the following result: 

Theorem 4.1 (Vanishing viscosity limit in phase dynamics). Suppose 
that the first order system (|1.2[) is uniformly hyperbolic for large value ofw 
hut may admit elliptic regions in the phase space, that is for some constants 
M, Co > 

inf o^v{w) > Cq-. 

\w\>M 

Suppose also that there exist constants Ci > and i] > such that \Ow{w)\ < 
Ci\w\'^~''' for \w\ > 1. Then for all Riemann data VyW; and Vr,Wr in the 
hyperbolic region of the phase space, i.e. 

Wi, iVr & "K := \w I Ow > o| 

the viscous Riemann problem (|1.4|) - (|1.5|) , with 6 = 0, admits a solution 



which has uniformly bounded total variation at least away from 
y = 0, and more precisely 

I \\v'^\ + y \w'J\\ dy < - Vi\ + \Wr - Wi\. 
Jr 

The functions Ve,We converge pointwise at all y Q to a limit v,w, which 
is a solution of the Riemann problem (|1.2|) - (|1.3|) away from the axis y 



xjt = Q. Furthermore, v has bounded variation and so admits left- and 
right-hand limits at y = 0, while the function w and its variation measure 
dw/dy satisfy 

, , 1 

lyl 

, , , , 
lyl li-l < oo. 

R dy 

Hence, the component v only has globally bounded variation. In 
turn, the conservation laws (|1.2[) are satisfied in the two regions x < 
and X > 0, but a singularity may arise on the axis. 

Proof. Step 1. A priori estimates for the components v, w. To simplify the 
notation we suppress the subscript e. We will first study the problem 
away from the axis y = 0. We consider the function 

«(y) = y - Mwiy)) - e)ly, (4.1) 
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and, given some 6 > 0, we study the problem in two regions [-L, -6] 
and [5, L], with the following boundary conditions 

w{-L) = ivi, w{6-) = w», (4.2) 

w{6) - w^, w{L) = Wy. (4.3) 

We set 

"P-^^^ J-5^-(mf;_aHx)dx' ^'^^^^ •= jL^-(l/e)f;;^a'Hx)dx' 

/y 
a{x)dx attains 

its global minimum. 

Let us consider A^"^ = sup^<^<^ {Ow{w) - e)"^^^^ and let c > max{l, A^"^}. 
For y > c 

{x^ - o^,{x) + e)dx > — -dx 

Jp+ Jc X 



)dx 



{x - c)dx 



2 

So, for y > (1 + sup^^ ^^((c7^(zi;) - e)*^'^, we get 



(piy)dy < -e-— (4.4) 



and ^ 

I y(piy)dy < (2 + sup ((ct^(w) - eyf\ (4.5) 



We then set 

w{y) = 



zvi + (w, - zvi) J^^ (pt dx, y < -6, 
Wr + (w, - Wr) (pi dx, y > 6. ^^'^^ 



Here, we have taken w(-6) = w(6) = w» and we also take v{-b) = 
v{6) = V*. Integrating the equation v' = yiv' from -L to -6 and 6 to L 
and using (|4.6|) , we get 

v^-vi = (w, - wi) j -y(p-{y)dy 

and 



- U» = (W» - Wr) 



y<p+iy)dy. 
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Adding these formulas, we get 

Vr-vi + iVr ycp^dy - wi y (p^_ dy 



w» := 



£y(pidy- ^''ycptdy 
and, subtracting, 

[Vr -vi + {wr - wi) 1/ (pi dy] -y (pt dy 



(4.7) 



= Vi + 



J^yfUy-f'^ycptdy 



(4.8) 



Now, V can be expressed in terms of the boundary data Wi, Wr, Vi, Vr 
and 6: 



^ [v,-vi+(wr-wi) J^yyj dy] ^\ -xipl dx ^ _^ 



Vr + 



[vi-Vr+(iVr-ivi) J_ L^-yipl dy] ^ x(p'\ dx 



(4.9) 



X yv'Uy+f^^ -y(p-dy 

Also, w{y) can be written as 

f'y.dx, 



L> y > 6. 



w\y) 



zvi + (z;» - vi) 



Wr + {Vr - 'V*) — 



ycpt dx 



y < -6, 



y>6. 



(4.10) 



Now clearly is bounded independent of e > 0,6 > : 

|U, - Vi\ < \Vr - Vi\ + L{\Wr - W/|) 

and from (|4.9|) and (|4.10|) , we get the following estimates 

\y\ \w\y)\ < {\vr\ + \vi\) + 2L(K| + \w,\), 
|z;^(y)| < 2{\vr\ + \vi\) + L{\wr\ + |w/|). 

Also we have the following estimates on the derivatives of v^, w^. 



X 



\yw\yy\dy < Qv,. - vi\) + LQw,. - W;|), 



v\yy\dy < {\vr - vi\) + L(K - W/|). 
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Furthermore, for all L > |i/| > 6 

\v\y)'\ < —[{\Vr - V^\) + L{\Wr - Wi\)y), 

\lv\y)'\<^J^{\Vr-Vi\) + L{\Wr-W,\)). 

The existence of solution for fixed L follows from these estimates as 
we will see. In order to pass L — ^ oo, we need estimates independent 
of L. For this we use the growth condition on Oy^. From the expression 
(|48)) for V,, and (|431) , we get 

\V* - Vl\ < \Vr - Vl\ + (2 + sup {Oioiw) - eY^^^)\Wr - Wl\. 

l<i/<L 

Using this in (l4!9b and RIOD , we get 

sup \w\y)\ < {\Vr\ + \vi\) + ((3 + sup (cj,„(w) - eY^'\\Wr\ + \Wi\), 

i<y<L 'i^iy^L 

\v\y)\ < 2{\vr\ + \vi\) + ((2 + sup {a,,{w) - e)+'^')(|zy,| + Ml 

l<y<L 

Now using our assumption |cri„| < Ci|wp~'' for |zy| > 1 in the first in- 
equality, we conclude that sup-^^^^^ l^t'^(y)l is independent of L. Sim- 
ilar estimate holds for [-L < y < -1]. We get there exists a constant 
C = C{vi, Vr, ivi, Wr), independent of L such that 

sup \w\y)\ < C (4.11) 

i<\y\<L 

With this constant C, let = sup\-^„\<c\Ow{w)^^^ . From (|4.4|) , we get 
(p^ is essentially supported in [6, -I- 1]. Similar arguments give cpt 
is essentially supported in the interval [-A^ - 1, -6]. It follows from 
(|4!9|) and (|4.10|) that there exists a constant Ci 



( |l/H;^(l/)'|rfy < {\Vr - Vi\) + Ci{\Wr - Wi\), 

J\y\>5 



j 

J\v\: 



(4.12) 



\v\yy\dy < {\vr - Vi\) + Ci(K - Wi\). 

'\y\>5 

Step 2. Existence proof. For each 6 > 0, we can apply Schauder 's fixed 
point theorem, as was explained earlier, and we obtain a solution 
(w^, u'^) defined on the interval [-L, -6] {J[5, L] and satisfying uniform 
in e > total variation estimate. The estimates (|4.11[) and (|4.12|) allows 
as to let L tend to infinity. By compactness, we have a solution {v, w) 
defined in the region \y\ > 6. Since 6 > is arbitrary, we obtain a 
well-defined solution away from y = 0. Indeed, has a uniform 
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total variation on R and v admits left- and right-limit at y = 0. One 
can check easily that the limit is a weak solution of the problem (|1.2[)- 
(|1.3|) away from the axis y = at least. Without a control on o{w) near 
y = we can not exclude a concentration term on the axis. □ 

5. Boundary Riemann problem and further generalizations 

In this section, we outline how the method developed in previous 
sections can be generalized to the boundary- value problem and in- 
dicate several generalizations. It is well-known that in the strictly 
hyperbolic case with nondegenerate function o, the boundary Rie- 
mann problem (in x > 0, f > 0) 

Vt - o{w)x = 0, Wf-Vx = 0. (5.1) 

with initial and boundary conditions 

w{0+,t) = Wh, v{x,0) = Vr,w{x,0) = Wr (5.2) 

is well posed. This easily follows from an analysis of the wave curves 
for the system corresponding to the left-moving and right-moving 
characteristic families. The physical regularizations considered in 
earlier sections is well-suited to handle this boundary value problem, 
without producing boundary layers in the limit. 

Consider the nonlinear elastodynamics with physical viscosity 

- yv' - o-,u{w)) w' = ev" , -yw' -v' = 0, on [0, oo), (5.3) 

with boundary conditions 

w(0) = Wh, v{oo) = Vr, w{oo) = W,: (5.4) 

Observe that, for the "full" viscosity approximation, we need to 
prescribe the component v at y = and this generates a boundary 
layer at y = 0, after passage to the limit e — > 0; see l|TOl [TT| . In the 
present case, we show that no boundary layer arises. 

As in Section 2, the problem can be reduced to a scalar equation 
for the unknown iv 

{y^ + e - Ozv{w)) w' = -ey w", on [0, L], (5.5) 

for sufficiently large L, with boundary conditions, 

w{0) = iVh,iv{L) = iVr. (5.6) 

Once we have w, we get the component v from the equation yiv' +v' = 
and the boundary condition v{L) = v,-. 

The fixed point argument in Section 2 yields a solution of (|5.5|) - 
(|5.6)) which is of uniformly bounded variation. Furthermore, it can 
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be represented by an integral formula in terms of the function <p+ 
introduced in Section 2 

iv'^iy) = Wy + {uo\, - w,-) I <p+ dx, y > 0. 
Using that (p+{x)dx = 1, we get 

r 

w'^iy) = Wb + {wh - Wr) I <p+ dx, y > 0. (5.7) 

Jo 

Using (|5.7|) in the equation v' = -yw' and integrating from y to L, we 
get ^ 

v^{y) = Vr + {wh - Wr) I x(p+ dx, y > 0, (5.8) 

where we used v{L) = v,-. Note that we cannot prescribe boundary 
condition at u(0) arbitrarily, since from the above equation it follows 
that ^ 

z;(0) = Vr + {wb - iVr) I x(p+ dx. 

Jo 

As in Section 2, using the fact that (p+ decays exponentially we can 
let L tend to oo. Again, using the properties (|2.9|) of cp+ near y = in 
(|5.7|) , it follows easily that there exists C > 0, a constant independent 
of e, such that for all y > close to the origin 

\w{y) -Wb\ <Cy, 

Hence, no boundary layer arises in this approximation at y = and in 
the limit as e — > 0, the limit function w = iv{y) satisfies the boundary 
condition w{0+) = Wb- We summarize our results in the following 
theorem: 

Theorem 5.1. Assume that infow > > 0. Then given arbitrary Rie- 
mann boundary data, the problem (|5.3|) - (|5.4|) admits a solution , which 
has uniformly bounded total variation and converges to a limit v, w which 
is a solution of the boundary Riemann problem (|5.1|) - (|5.2[) . 



We now discuss a general system. The basic nature of the system 
(|5.1|) is that V appear linearly and the characteristic speeds are equal in 
magnitude and opposite in sign. Further using the second equation, 
the system can be reduced to uncoupled equations for w. After 
this reduction we solved for w first and then for v. Our results 
can be generalized for systems of first order equations for w and v 
vector valued functions, having same structure, using the ideas of 
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the present work and the earlier work |[TOl ITT]] on boundary value 
problems. 

For example let F : — > be a smooth function with A{w) = 
D^F has real distinct positive eigenvalues Ai(w)^ < A2(w)^ < ... < 
Ajv(w)^ with a complete set of right-left normalized eigenvectors rj, j = 
1,2,. ..N, Ipi = 1,2,. ..,N, lj.n = 6jk. 

We assume that there exists Cq > such that Aj{w) > Cq for all 
j = 1,2, ...N and for all w e B{5o), a ball of radius 6o around a fixed 
state that we take to be 0. We consider the system of 2N equations 

(5.9) 

aw _ 

where v = v{t, x) and iv = iv{t, x) are valued functions. The system 
(|5.1|) is a special case of (|5.9|) with N = 1 and A(w)^ = c7t,;(w) 

We consider boundary value problem for (|5.9|) on x > 0, i > 0. 
When physical viscosity terms are taken into account, the corre- 
sponding set of differential equations becomes 

-yv'-F{wy =ev", 

- yw - V = 0, 

supplemented with the boundary conditions 

\im{v, w){y) = {Vr,Wr), 'w{0+) = Wh- (5.11) 

We first consider (|5.10|) on [0, L] with boundary conditions 

{V,IV){L) = {Vr,lVr), 'W{0+) = Wh- (5.12) 

As in the previous case the problem can be reduced to first solving a 
system for w{y) namely 

iy^ + e- A{iv)) w' = -ey w" , (5.13) 

on [0, L], with boundary conditions, 

w(0) = iVh,w{L) = Wr. (5.14) 

Existence of uniformly BV solutions of this problem easily follows 
from the work of Joseph and LeFloch IITOl [TT| . We just outline the 
main steps omitting the details. 
We decompose u' in the basis of eigenvectors of A{u), 

N 

u'iy) = '^^ajrj{u),aj = {lj,u') (5.15) 

7=1 
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A straight forward calculation lead to the following system of non- 
liear equations for Uj, j = 1,2, ...,N 

eya'j + {y^ + e- A^{w))aj = Di{a,a) (5.16) 

where 

Di{a, a) = -ey aiUkiDuriSk) (5.17) 

k,i 

Observe that the linearized form of the equation (|5.16|) has the form 

e y a'jiy^ + e - Aj{w)) a'j = 0, 

and the corresponding wave measure (pj+ has exactly same proper- 
ties (|2.9|) given in section 2. These wave measures naturally appear 
when we invert the linear part of the equation (|5.16|) . Because of the 



quadratic righthand side this leads to interaction terms with differ- 
ent families of wave measures which is controlled by Tj%i (Pj+- Fixed 
point arguments give a uniform BV solution w^iy) for the system 
(|5.13)) and (|5.14[) . The details are similar to Joseph and LeFloch IITOl 



and, therefore, are omitted. Then the BV estimate for v follows from 
the second equation in (|5.10|) and 



v{y) = Vr+ I xw'{x)dx. 

Jy 

As the estmate (|2.9|) shows that (pj+ is essentially supported in 
[A J', Aj^] we can let L go to infinity. Here AJ' and A^ denotes the 
minimum and maximum of Aj{w) on the ball B(6o). We have the 
following result. 

Theorem 5.2. Assume that inf Aj{w) > Cq > 0. Then there exists bi > 
0,62 > such that, for Wh,iOr e B(6i) and arbitrary Vr the boundary 
Riemann problem (|5.10|) - (|5.11|) admits a solution , which has uniformly 



bounded total variation with w^{y) e 6(62) and converges to a limit v,w 
which is a solution of the equation (|5.9|) with boundary conditions w{0+) = 

Wh, W{X, 0) = Wr, V{X, 0) = Vr- 

Remark 5.3. This analysis is easily extended to the case with physical 
viscosity and capillarity as well namely 

-yv' - F{wy = ev" - ye^w'" , y > 0, 

, \\ ^ ^ (5.18) 

-yw -V =0, 

supplemented with the boundary conditions 

\m\{v,w){y) = {Vr,w,), w{0+) = W\,. (5.19) 
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As in the previous case, the problem can be reduced to first solving 
a system for w{y), namely 

iy^ + e- A{w)) w' = -ey w" - ye^w'" , (5.20) 

on [0, L], for sufficiently large L, with boundary conditions 

w(0) = iVh,iv{L) = iVr. (5.21) 

We decompose u' in the basis of eigenvectors of A{u), 

N 

"'(3/) = i^j^jiu), Uj = {Ij, u') (5.22) 

and then the nonlinear equations for Uj, j = 1, 2, N take the form 
ye^a'l + eya'- + {y^ + e- \^-{w))aj = Di{a, a) + D2{a, a) + Do,{a, a, a), (5.23) 
where 

Di{a,a) = -ey^aiUkiDuTi.rk) 

k,i 



D2{a,a) = -ye\^aia'f,{Duri.rk) + ^{aiakYiDuri.Tk 

k,i k,i 

D3{a,a,a) = -ye^^aiakaiDu{Duri.n)ri 



k,i,l 

The linearized equation for (|5.23|) is 

ye^aj + eya'- + {y^ + e- A^-{w))aj = 0, 

and hence the wave measures (pj+ in this case has the same qualitative 
properties as (|3.15|) . 

Existence of uniformly BV solutions of the problem (|5.20[) and (|5.21|) 
can be easily deduced from [17J. Then, the BV estimate for v follows 
from the second equation of (|5.18|) . The details are omitted. We get 
the following result: 

Theorem 5.4. Assume that inf Aj{w) > Cq > 0. Then there exists 6i and 62 
such that for Wb,Wr e B(6i) and for arbitrary Vr, the problem (|5.18|l - (|5.19|) 
admits a solution , , which has uniformly bounded total variation with 
w^iy) e 6(62) and converges to a limit v, w which is a solution of the equation 



(|5.9|) with boundary conditions w{0+) = Wb, w{x, 0) = Wr, v{x, 0) = v,-. 

Now let us consider the phase transition case with physical viscos- 
ity and capillarity. Following previous sections , the system can be 
reduced to a single scalar equation for the unknown w, 

(y^ + e - o,u{w)) w' = -ey w" - ye^w'" , 7 > (5.24) 
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which can be studied in two regions [-L, -60] and [-60, L] away from 
y = 0, with the boundary conditions 

w{-L) = Wi, w(-6o) = w», (5.25) 

w(6o) = w*, w(L) = Wr, (5.26) 

where w, is given by (|4.7|) with 6 replaced by 60 chosen as follows. 
Following Sections 3 and 4, we can construct a BV solution w for 
(|5.24|) , in [-L, -60] and [60, L] with the boundary conditions (|5.25|) 
(|5.26|) , provided 

1 e 
iU(w, y) = cT,„ + i/^(— - 1) - - > 0. 

This is the case if we impose the restriction \y\ > 60 where 

Here, c > is a constant such that Oy,{w) > -c for all w. 
We get the following result : 

Theorem 5.5 (Vanishing viscosity-capillarity limit in phase dynam- 
ics). Suppose that the first-order system (|1.2|) is uniformly hyperbolic for 
large value ofw hut may admit elliptic regions in the phase space, that is for 
some constants M,Cq> 

inf Ow{w) > Co^. 

\zd\>M 

Suppose also that there exist positive constants Ci,ri such that \a^\ < Ci|zfp"'' 
for |w| > 1. Let c > be such that o^{w) > -c for all w and suppose 60 
is chosen as in (|5.27[) . Then, for all Riemann data Vy Wi and Vr, iVr in the 
hyperbolic region of the phase space, i.e. 

zvi, zvr G "K := / Oi„ > o| 

the viscous-capillarity Riemann problem (|1.4[) - (|1.5|) , with 6 = ye^, y > 
admits a solution v'^,w^ which has uniformly bounded total variation in 
\y\ > 60, and more precisely 



I 

J\v' 



(\v'^\ + \w'^\j dy < \Vr - Vi\ + \Wr - Wi\. 

y\>5o 

The functions Ve, We converge pointwise at all \y\ > 60 to a limit v, w, which 
is a solution of the Riemann problem (|1.2 |) - (|1. 3)1 /or \y\ = \x/t\ > Sq. 
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For instance, the example o{w) = u{u^ - 1) satisfies all the as- 
sumptions of the theorem. We can not exclude that the vanishing 
viscosity-capillarity approximations to the Riemann problem could 
contain a highly oscillating (i.e. weakly but not strongly converg- 
ing) pattern near the axis; its thickness would be 26o, at most. The 
existence of such stationary waves is not surprising as it has been 
observed numerically; they were never pointed out analytically until 
now. 



6. Effect of general viscosity in fluid dynamics 

Following our earlier work [12J we now return to the general dif- 
fusion approximation for strictly hyperbolic systems with general 
diffusion matrix B{u) of the form 

- yu' + A{u)u' = e{B{u)u'y, (6.1) 

and we consider an associated generalized eigenvalue problem. 

Let /.ty, Ij, Tj, i = 1,2, ...,N be the eigenvalues and left and right- 
eigenvectors given by 

(-y + y) = l^jiv, y) B{v)7j{v, y), 

^ / . ^ (6.2) 

Let us impose the normalization 

li{v, y) B{v)7j{v, y) = if z ^ ;; lj{v, y) B{v)Tj{v, y) = 1. 
In the special case B{u) = I we have simply 

y) = -y + A;(u), 'Tjiv, y) = Tjiv), Ijiv, y) = Ijiv). (6.3) 

So, we expect the eigenvalues and eigenvectors in (|6.2|) to be close 
to to that of (|6.3|) , at least in the case that B{v) = I + rj T{v, rf) where 

T{v, 7]) = (tijiv, rj)) with tijiv, 7]) = 0(1) and t] « 1. 

Lemma 6.1. If \B{v) - 1\ = 0(7]) is sufficiently small, then 

^i{'v,y) = -y + A;(u) + 0(ij), lj{v,y) = Ijiv) + 0(r]), 

_ - _ (6.4) 

Tjiv, y) = Tjiv) + Oii]), li{v, y) B{v) dyVjiv, y) = 0(r]). 

Proof. When B has the form B{v) =1 + 7] T{v) where the matrix T{v) 
does not depend on rj, the desired estimates follow from classical 
results. In the slightly more general case when T = T{v, rj) we may 
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argue as follows. First, note that the eigenvalues under consideration 
are given by 

det[[-yI + A{v))-^B{v)) = 0. (6.5) 

Clearly, (|6.5|) is a polynomial equation in [j. whose coefficients are 
polynomials in the i/y's. The leading coefficient det(l + i] {tij^^ is dif- 
ferent from zero. Since the roots of (|6.5|l are distinct for tij = 0, then 
for i] sufficiently small the roots of (|6.5|) (namely the /.zy's) are distinct 
and depend smoothly upon the 's with, in addition, 

^ji^r y) = -y + A;(u) + riYj ^iji'"' y) + (6-6) 

i 

where the coefficients are smooth functions in y. 

Since the corresponding left- and right-eigenvectors (which need 

not be normalized at this stage) Ij and Tj are polynomials in and 
it follows that they are also smooth in the fl,/s. Substituting them 
in (|6.3)) we get 

hi'^' y) = hi'^) + ^'i^'^' ']) + •••' 

^ (6-7) 
Tjiv, y) = Tjiv) + T]2_^ Tijiv, y) tij{v, 1]) + ..., 

t 

where Uj and r;y depend smoothly upon y. From (|6.7|) , it follows 
that for rj sufficiently small the expansion holds for the normalized 



vectors lj{v, y) and rj{v, y). Finally, the required estimates (|6.4[) follow 



from (|6.6|) and (|6.7|) . This completes the proof of (|6.4|) . □ 

By the same technique as above we can also prove: 

Theorem 6.2. Under the conditions that \DuB{u)\ < r] and the characteris- 
tic fields associated with A{u) are genuinely nonlinear, the generalized Lax 
shock inequalities 

Xjiu{y+),y) <y< Aj{u{y-),y) 
are equivalent to the standard Lax shock inequalities, 

Aj{u{y+)) <y< Aj{u{y-)), 

where Ay = (xj,ATj). 

Proof. We assume that VAy • ry > 0. The key observation is that 

VAy • Tj is also positive and therefore the same part of the Hugoniot 
curve is selected by the standard and by the generalized Lax shock 



PHASE DYNAMICS WITH PHYSICAL VISCOSITY AND CAPILLARITY 27 

inequalities. This is so because VAy ■ Tj = VAj ■ tj + 0(r]) as we now 
check. 

When |D„B(i/)| < r]weseethatDi,r(z;, < C and that the expansion 
\6.7) is valid for the u-derivatives as well. This completes the proof. 

□ 

We end this section with a further discussion of systems of two 
equations, and provide some explicit calculations which lead to suf- 
ficient conditions on the diffusion matrix B allowing us to apply the 
techniques introduced earlier in this paper. Note first that (|6.5)) can 
be written as 

det((-yJ + diag(Ai(z;),A2(u)) - ^L{v)B{v)R{v)) = 0, (6.8) 

where L{u) and R{u) are matrices of left- and right-eigenvectors as- 
sociated with A{u). Namely, L{v) := {li{v), . . . ,l^{v)^ and R{v) := 

(ri{v), TNiv)). Setting L{v) B{v) R{v) = becomes 

(detB(i;)) ^2 - (bniv) {A2 -y) + b22{v) (Ai - y)) fi + (Ai - y) (A2 - y) = 0. 
Solving this quadratic equation in [j. and using the notation 

. biijv) b2i{v) , , Hv) - y 

we arrive at 

(a2{v, y) + a,{v, y)) - (a2{v, y) - aM y)) (l + 4 f^iv) g^H)''' 
^'^""'^^ = 2(1^) ' 

(a2{v, y) + aM y)) + (a2{v, y) - a,{v, y)) (l + 4^(z;) ^T^ '^ 

^'^''^^ = 2(1^) • 

It is easy to see that and ^2{v) are real if < (i{v) < 1, which is 
guaranteed if 

bn{v) > 0, b22{v) > 0, 

bi2{v)b2i{v)>0. ^^-^^ 

The first condition above is also a necessary condition for viscous 
shocks to be strictly stable in the sense of Majda and Pego 1120 1 . 

Assuming (|6.9|) , an easy calculation shows that the eigenvalues 
jUi(u) are separated: 

M^) - M^) > C > 0. 
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It is easy to compute the corresponding right- and left-eigenvectors, 
namely 

Ti{v, y) = [- bi2{v) ^i{v, y), Ai{v) - y - hii{v) ^i{v, y)) , 
k{v, y) = ( - &2i(u) iii{v, y), Ai{v) - y - hii{v) ^i{v, y)). 



Finally, it is easily seen that the estimates (|6.4|) hold true when |B(w) - 

/| = o{^). 
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